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ABSTRACT 

Context. Pure gas-phase chemistry models do not succeed in reproducing the measured abundances of small hydro- 
carbons in the interstellar medium. Information on key gas-phase progenitors of these molecules sheds light on this 
problem. 

Aims. We aim to constrain the chemical content of the Horsehead mane with a millimeter unbiased line survey at two 
positions, namely the photo-dissociation region (PDR) and the nearby shielded core. This project revealed a consistent 
set of eight unidentified lines toward the PDR position. We associate them to the 1-C3H + hydrocarbon cation, which 
enables us to constrain the chemistry of small hydrocarbons. We observed the lowest detectable J line in the millime- 
ter domain along a cut toward the illuminating direction to constrain the spatial distribution of the 1-C3H + emission 
perpendicular to the photo-dissociation front. 

Methods. We simultaneously fit 1) the rotational and centrifugal distortion constants of a linear rotor, and 2) the 
Gaussian line shapes located at the eight predicted frequencies. A rotational diagram is then used to infer the exci- 
tation temperature and the column density. We finally compare the abundance to the results of the Meudon PDR 
photochemical model. 

Results. Six out of the eight unidentified lines observable in the millimeter bands are detected with a signal-to-noise 
ratio from 6 to 19 toward the Horsehead PDR, while the two last ones are tentatively detected. Mostly noise appears at 
the same frequency toward the dense core, located less than 40" away. Moreover, the spatial distribution of the species 
integrated emission has a shape similar to radical species such as HCO, and small hydrocarbons such as C2H, which 
show enhanced abundances toward the PDR. The observed lines can be accurately fitted with a linear rotor model, 
implying a E ground electronic state. The deduced rotational constant value is B — 11244.9512 ± 0.0015 MHz, close 
to that of I-C3H. 

Conclusions. This is the first detection of the 1-C3H + hydrocarbon in the interstellar medium. Laboratory spectroscopy 
is underway to confirm these results. Interferometric imaging is needed to firmly constrain the small hydrocarbon 
chemistry in the Horsehead. 

Key words. Astrochemistry - ISM: clouds - ISM: molecules - ISM: individual objects: }{ orse }i ea( i nebula - Radio lines: 
ISM 



1. Introduction 

Molecular ions play an important role in the physics and 
chemistry of the insterstellar medium. They trace the gas 
physical conditions, e.g., its ionization rate and its ioniza- 
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tion fraction (see for example iGoicoechea et all 12009). and 
they participate in the coupling of the gas with the mag- 
netic field. Moreover, molecular ions are key species in the 
gas phase synthesis of molecules because ion-molecule re- 
actions most often have no activation barrier. 

Simple hydrocarbon molecules, such as C2H, C3H, and 
C3H2 have been detected in a wid e variety of sources 
from diffuse (e.g. iLucas fc L iszt 2000) to dark c louds (e.g. 



IWootten et al"1ll980t iManeum k Woottenlll990D . The high 
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Fig. 1. Millimeter lines attributed to the I-C3H" 1 " cation in the PDR (lower panels) and the dense core (upper panels) 
positions. The numbers above each panel are the line rest frequencies in MHz. The spectra and the Gaussian fits are 
shown as black histograms and green curves, respectively. The inferred rest frequencies are displayed as vertical dotted 
black lines for v\ sr = 10.7 km s -1 and ±3ct noise levels as horizontal dotted red lines. 



abundances found at the UV-illuminated edges of molecu- 
lar clouds or photo-dissociation regions (PD Rs) cannot be 
reproduced by curren t pure g as-phase mo d els (iFuente et al.l 
I2003t iTevssier et all I2004D . iPetv et alj ((20051) proposed 
another chemical route, namely the photo-erosion of 
Polycyclic Aromatic Hydrocarbons (PAHs) and small car- 
bon g rains to produce these hydrocarbons. iRimmer et al.l 
(|2012h showed that using a column-dependent cosmic 
ray ionization rate slightly improves the agreement be- 
tween models and observations in the Horsehead PDR. 
Ion-molecule reactions with the C2H 4 " and CsH + cations 
are thought to be the most import ant gas-phase chan- 
nels to form smal l hydrocarbons (jTurner et al.l l2000t 
Wakelam et al.l[2010h . but they have not yet been observed 
in the interstellar medium. Constraining the abundances of 
these intermediate hydrocarbon cations will shed light on 
the formation routes of hydrocarbons. 

In this paper, we report the first detection in the inter- 
stellar medium of the I-C3H" 1 " hydrocarbon cation toward 
the Horsehead PDR (more prec isely, at the peak of the 
HCO emission, Gerin et al.ll2009D . Sect. [2] describes the ob- 
servations. Sect. |3] and @] explain how we inferred the spec- 
troscopic parameters associated with the set of unidentified 
lines, and why we attribute these to 1-C3H + . Sect. [5] dis- 
cusses the determination of the 1-C3H + abundance and its 
chemistry. 

2. Observations 

The Horsehead WHISPER project (Wideband High- 
resolution Iram-30m Surveys at two Positions with Emir 
Receivers, PI: J. Pety) is a complete unbiased line survey 
of the 3,2, and 1 mm bands that is being completed at the 
IRAM-30m telescope. Two positions are observed: 1) the 
"HCO peak" (RA=5h40m53.936s, Dec=-2°28'00", J2000), 
which is characteristic of the phot o-dissociation regi on at 
the surface of the Horsehead nebula (|Gerin et al.ll2009l ). and 
2) the "DCO+ peak" (RA=5h40m55.61s, Dec=-2°27'38", 
J2000), which belongs to a cold and shielded condensa- 
tion located less than 40" a way from the PDR edge, where 
HCO+ is highly deuterated (|Petv et al.ll2007h . The combi- 
nation of the new EMIR receivers and Fourier transform 
spectrometers at the IRAM-30m telescope yields a spec- 
tral survey with unprecedented combination of bandwidth 



(36 GHz at 3 mm, 19 GHz at 2 mm, and 76 GHz at 1mm), 
spectral resolution (49 kHz at 3 and 2 mm; and 195 kHz at 
1mm), and sensitivity (median noise 8.1 mK, 18.5 mK, and 
8.3 mK, respectively). A detailed presentation of the ob- 
serving strategy and data reduction process will be given in 
another paper. In short, any frequency was observed with 
two different frequency tunings and the Horsehead PDR 
and dense core positions were alternatively observed every 
15 minutes. The total observing time amounted to one hour 
per frequency setup and position. 

In the analysis of the survey, we found an 89 mK line 
peak around 89.957 GHz, that could not be associated 
to any t ransition listed in the com mon public catalogs : 
CDM$| (iMuller et alJ l2001h . JPlQ i|Pickett et all Il998h . 
and splatalogufl The observing strategy allowed us to rule 
out that the detected line is a ghost line incompletely re- 
jected from a strong line in the image side band (the typical 
rejection of the EMIR sideband separating mixers is 13 dB). 
Our search for an identification started with the simplest 
assumption, i.e., the associated species is a linear rigid ro- 
tor. The frequency of the transition ( J + 1 J) is then 
given by 

u = 2B(J + l), (1) 

where B is the rotational constant. In this simple model, 
the ratio of line frequencies depends only on the ratio of 
J + 1 values. We thus could predict different sets of fre- 
quencies to search for companion lines associated to the 
same species, each set associating a given J + 1 — > J transi- 
tion to the frequency of the detected unidentified line, i.e., 
~ 89.957 GHz. The only set of frequencies that consistently 
brings five other detected lines and two more tentative de- 
tections less than 1 MHz from the frequency predictions was 
the one which associates the 89.957 GHz unidentified line to 
the J = 4^3 tr ansition of a linear rigid rotor. Using the 
WEEDS ext ension fMaret et al.l 1201 ll ) of the GILDAS/CLASS 
softwar^l (Pety 2005), we quickly ruled out the assignment 
of any of the detected unidentified lines to other possible 



http : / /www . astro . uni-koeln . de/ cdms/ 

2 http://spec.jpl.nasa.gov/ 

3 http://www.splatalogue.net 

4 See http://www.iram.fr/IRAMFR/GILDAS for more informa- 
tion about the GILDAS softwares. 
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Table 1. Results for the simultaneous fit of the 1-C3H+ line parameters (bottom) and the associated spectroscopic 
modeling (top). 
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a C is the fit likelihood. b Parameters for a linear rotor developed to second and third orders (see Sect. [2]). c Results of the 
simultaneous Gaussian fits using the third-order model to predict the frequencies. The common line width is 0.81 ± 0.03 km s~ 

d Outside of the EMIR receiver tuning range. 



species because the potential candidate species were com- 
plex molecules for which many other expected transitions 
were not detected in the survey. 

Fig. [1] displays the spectra at the Horsehead PDR and 
dense core positions of the consistent set of eight uniden- 
tified lines, which lie in the millimeter frequency bands. 
Even though the weather conditions, the pointing and fo- 
cus corrections, and the tuning setups were shared for the 
two observed positions, the unidentified lines were detected 
with a profile, integrated signal-to-noise ratio between 3 
and 19 at the PDR position, while only two-sigma upper 
limits of typically 25mKkms" 1 for a 0.8 km s -1 linewidth 
could be derived at the dense core position. The only ex- 
ception is the J = 4 — 3 line, which is tentatively detected 
toward the dense core. Assuming that the emission arises 
from a Gaussian filament of 12 " full width at ha lf maxi- 
mum centered on the PDR fsee lGerin et al.ll2009h . 27% of 
the emission detected at the core position is explained by 
beam pickup from the PDR. The remaining emission could 
arise in the low er density skin of t he dense c ore, already de- 
tected in HCO (jGerin et al.ll2009h and CF+ (|Guzman et al.l 
l2012allbT ). In summary, the eight (tentatively) detected lines 
in the PDR are unlikely to be observing artifacts. 

3. Associated spectroscopic constants 

To compute the spectroscopic parameters associated with 
the set of unidentified lines, we used higher order improve- 
ments to the simple model of a linear rotor. The second 
and third order corrections, which include the effect of 
the centrifugal distortion, predict that the frequency of the 
J + 1 — > J transition is 

v = 2B(J+1) -4L>(J±1) 3 , and (2) 

v = 2B(J+l) — 4:D(J+lf + H(J+ 1) 3 [(J±2) 3 — J 3 ], (3) 

where D and H are the centrifugal distortion constants to 
the second and third order, respectively. 

Using these expressions, it is possible to generate the 
sum of 1) Gaussian white noise, 2) a baseline offset, and 



3) Gaussian line profiles at the eight predicted frequencies 
of arbitrary areas and line widths. This allows us to simul- 
taneously fit the noise level, a residual baseline offset, and 
the Gaussian parameters (areas and line widths) as well as 
the rotational and centrifugal distortion constants associ- 
ated with these models. In practice, it is very difficult to 
make the fit converge with the eight lines because of the 
limited signal-to-noise ratio of some of them. We therefore 
assume that the eight lines have the same line width. This 
is likely for two reasons. First, the line width is dominated 
by the turbulent velocity field and the lines are probably 
emitted from the same gas cells because the range of energy 
probed by the eight transitions is relatively narrow (from 10 
to 80 K). Second, the velocity gradient in the plane of the 
sky is shallow, ensuring that the line width does not evolve 
significantly when the beam size changes with frequency. 
This assumption enables a quick fit convergence when we 
start the fit with the following initial values: B = 11.2 GHz, 
D = 5 kHz, H — 1 Hz, zero offsets, 10 mK noise levels, and 
1 kms" 1 common linewidth. Appendix [X] discusses why the 
third order model gives the best fit to the data and Table [T] 
gives the results of the fits (both the Gaussian line param- 
eters and the spectroscopic constants). 

Because laboratory measurements are not available, the 
only unambiguous frequency available is the one measured 
in the local standard of rest (LSR) frame, v\ sr . The fre- 
quency given in the source frame, ^ sou , assumes an LSR 
systemic velocity of the source, Ui sx . These frequencies are 
linked through i/ sou = v\ sv [1 — (v\ sr /c)] , where c is the 
speed of light. The systemic velocity may vary from species 
to species because of different coupling between the gas 
kinematics and its chemistry and/or line-shape variations 
due to optical-depth effects. The analysis of the unambigu- 
ously attributed detected lines in the survey shows that the 
possible range of LSR systemic velocities in the Horsehead 
PDR is quite narrow, ~ 0.2 kms -1 around 10.7 km s^ 1 (see 
Fig. IA.1I which displays different lines of several hydrocar- 
bons). All frequencies quoted here accordingly assume an 
LSR velocity of 10.7kms _1 . The derived spectroscopic pa- 
rameters may therefore need to be slightly linearly scaled. 
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Fig. 2. High angular resolution maps of the integrated intensity of H 13 CO + , HCO, C2H, and C-C3H2, and the 7.7 /jm 
PAH emission. Maps are rotated by 14° counter-clockwise around the projection center, located at (5x, Sy) = (20", 0"), to 
bring the illuminated star direction in the horizontal direction. The horizontal zero is set at the PDR edge. The emission 
of all lines is integrated between 10.1 and 11.1 km s -1 . Displayed integrated intensities are expressed in the main beam 
temperature scale. Contour levels are shown in the grey-scale lookup tables. The red vertical line shows the PDR edge 
and the blue and green crosses show the dense core and PDR positions, respectively. 
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Table 2. Lines simultaneously observed in the cut dis- 
played in Fig. [21 
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Fig. 3. Integrated intensity profiles across the PDR photo- 
dissociation front of the I-C3H+, HCO, H 13 CO+, and C 2 H 
species (IRAM-30m). The cut would appear horizontal at 
Sy = 0" on Fig. [5] The observed lines are summarized in 
Table H 



4. Attribution of the lines to the l-CsH + cation 

The unidentified lines were all detected at the PDR po- 
sition, but not at the dense core position. We thus com- 
plemented these data with a cut from the Hn region into 
the molecular cloud along the direction of the exciting star, 
a Ori, i.e., perpendicular to the photo-dissociation front. 
Figure [3] displays the integrated intensity emission of the 
J = 4 — 3 unidentified line as a function of the angular 
distance from the photo-dissociation front, along with the 
emission of HCO, H 13 CO+, and C 2 H species (Tabled de- 
fines the observed lines). All these lines were observed si- 
multaneously with the IRAM-30m telescope during 4 hours 
of mild summer weather (typically 11mm of precipitable 
water vapor), but they all were well-detected nevertheless. 
At 29" resolution, it is clear that the unidentified species cut 
peaks in the UV-illuminated part of the Horsehead mane 
in the same way as HCO and C2H, while H 13 CO + peaks 
in the dense core, shielded from the UV field. This led us 
to conclude that the species is a reactive molecule with 
a spatial distribution similar to small hydrocarbon chains. 
For reference, Fig. [2] displays the integrated intensity emis- 
sion of the same lines, imaged at ~ 6" with the Plateau 
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H i3 CO+ 
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86754.2884 
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86777.4600 


C 2 H 


N = 1 - 0, J = 3/2 - 1/2, F = 2 - 1 


87316.8980 


I-C3H+ 


J = 4-3 


89957.6250 



de Bure Interferometer (jPetv et al.ll2007t iGerin et al.ll2009fc 
iPetv et al.| [ 2005t). plus the 7.7 umPAH emission imaged at 
~ 6" with ISO (Abergel et al.ll2003r) . 

The quality of the spectroscopic fit suggests that the 
molecule is a linear rotor with a 1 S + electronic ground 
state, i.e., with a closed electronic shell. The rotational 
constant is B ~ 11.24 GHz, which implies the presence of 
several heavy atoms in the species. According to the lit- 
erature, the most probable canditate is the 1-C3H + cation. 
Indeed, ab initio calculation implies that 1) the linear struc- 
ture is the most stable, 2) it has the right electronic state, 
and 3) th e computed rotational constant value is about 
11.1 GHz dRadom et al.lfl976L IWilson fc Greed [l980L [19821: 



ICoooer fc Murphvlll98a Ilkutalll997h . Experimental spec- 
troscopic confirmation for this cation is being performed at 
the PhLAM laboratory in Lille (Bailleux & Margules, priv. 

co mm.). 

IWilson fc Greenl (|1980l ) estimated the dipole moment of 
this cation to be 2.6 Debye, but this is a rather old value. We 
therefore computed it again with more sophisticated ab ini- 
tio techniques of quantum chemistry, as implemented in the 
MOLPRO suite of programfS We used the CASSCF-MRCI 
level of theo ry with the correla t ion-co nsistent aug-cc-pVQZ 
basis sets of lWoon fc: Dunni ng (1993) for all atoms. The ac- 
tive space of the CASSCF included the n = 2 orbitals of 
carbon and the IS orbital of H. The geometry was first 
optimized at the CASSCF level, leading to bond distances 
for I-C3H+ (1.090A, 1.2 46 A, 1 . 355 A ) that perfectly match 
the results obtained in llkutal (jl997h . The dipole moment 
was then computed at the center of mass of the molecule 
for the CASSCF optimized geometry. The resulting dipole 
moment of I-C3H+ is 3 Debye at the CAS SCF-MRCI level 
of theory, i.e., close to the initial value of Wilson fc: Greenl 



5 H.-J. Werner & P. J. Knowles, MOLPRO (version 2002.6) 
package of ab initio 254 programs, 2002. 
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Table 3. Comparison of the measured and modeled abundances with respect to the number of protons at three different 
positions across the PDR front. 
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a Petv ct al. (2005), deduced from observations at ~ 6"-resolution. 
b This work, using a simple source model to compute the beam dilution. 



T rot = 14±2 K 

N(m) = (mo/m) 2 4.8±0.9 10 11 cm" 
/2 = 3 Debye 



Table 4. Column densities and abundances with respect 
to the number of protons, i.e., [X] = 0.5 N(X)/N(R 2 ), 
toward the PDR position, from single-dish observations at 
resolutions between 25 and 28". 
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Fig. 4. Rotational diagram of 1-C3H+ at the PDR position. 

(|l980f ). We used the new value to compute the column den- 
sity of I-C3H 4 -. 

5. I-C3H 4 " abundance and chemistry 

We used the formalism of the rotati onal diagram to esti- 
mate the column density of 1-C3H + . Goldsmith & Langer 
(1999) give the expressions for the level energies, the quan- 
tum level degeneracies, the Einstein coefficients, and the 
partition function Z(T) for a rigid rotor of dipole moment, 
\i. For instance, the Einstein coefficients are given by 

. 647T 4 i/V J±l 

Aj+M - Shot 2J±3' (4) 

where h and c are the Planck constant and the light speed, 
respectively. Using these formula, we are able to compute 
the rotational diagram, assuming 1) that the lines are op- 
tically thin (r < 1), and 2) that the 1-C3H+ emission is 
co-spatial with the illuminated filament of the C2H emis- 
sion, i.e., it approximately fills as a Gaussian filament of 
~ 12" width in the 5x direction, infinite size in the 8y di- 
rection, and centered at the HCO peak (see Fig. [2]). The 
filling factors are 0.4, 0.6, and 0.8 at 90, 157, and 270 GHz, 
respectively. Figure0]displays the diagram for ^0 = 3 Debye 
(this work). The rotational temperature inferred from the 
observations is independent of the dipole moment because 
it scales with the slope of the rotational diagram, while 
the column density scales as fi~ 2 . The typical rotational 
temperature and column density thus are 14 ± 2 K and 
(Ho/fi) 2 (4.8 ± 0.9) x 10 n cm- 2 , respectively. 

We took advantage of the Horsehead WHISPER survey 
to make a consistent summary of the small hydrocarbon de- 
tections in the Horsehead PDR. The spectra of the detected 



lines in the 3mm band of C2H, C-C3H, I-C3H, C-C3H2, 1- 
C3H2, and C4H are displayed in Fig. I A. H and their Gaussian 
fit results are summarized in Table. IA.1I (both available on- 
line only). The integrated intensities of the lines are con- 
sisten t with the results published bv iTevssier et al. (2004, 
2005). We accordingly just su mmarized their colu mn densi- 
ties and abundances in Table l4"l [Petv et al~l(|2005l) presented 
higher angular (typically 6") observations of C2H and c- 
C3H2 obtained with the Plateau de Bure Interferometer. 
The associated abundances at three different positions (the 
"IR peak" close to the usual "HCO peak", one position 
named "cloud" representative of the UV-shielded material, 
and one position named "IR edge" closer to the Hn region) 
are summarized in Tabled The comparison of Table[3]andE] 
indicates that the C2H and C-C3H2 abundances measured 
at a typical resolution of 25 — 28" falls in the uncertainty 
range of abundances measured at a four times better reso- 
lution with PdBI. This means that their emission more or 
less fills the 25 — 28" beam, even though it is substructured 
into filaments (see Fig. [2). Indeed, if we had assumed that 
the 1-C3H + emission fills the beam at every measured fre- 
quency, the inferred I-C3H 4 " abundance would have only be 
reduced by a factor 2. 

We finally compared the derived abundances with a one- 
dimensional, steady-state photochemical model (Meudo n 
PDR model, iLe Bourlot et al.ll2012t iLe Petit et al1l2006l ). 
The used version of the Meudon PDR code includes the 
Langmuir-Hinshelwood and Eley-Rideal mechanisms to de- 
scrib e the formation of H2 on grains (|Le Bourlot et all 
120121 ). and surface reactions for other species (Le Bourlot 
et al., in prep.). The physical conditions in the Horsehead 
have already been constrained by our previous observa- 
tional studies and we kept the same assumptions for the 
steep den sity gradient, radiat ion field (x — 60 in Draine 
units, see lAbergel et al.l l20 03 f). elemental gas- phase abun- 
dances (see Table 6 in iGoicoechea et al.ll2006l) and cosmic 
ray prim ary ionization rate (C = 5 x 10 _17 s _1 per H2 
molecule, IGoicoechea et alll2009f ). 
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Fig. 5. Photochemical model of the Horsehead PDR displayed at the full resolution of the model (left column) and 
convolved with a Gaussian of 6"-FWHM (right column). Top: Horsehead density profile nn = n(H) + 2n(H2). Middle: 
Predicted abundance of selected small hydrocarbons molecules and cations. Bottom: Abundances of the same hydro- 
carbons relative to the abundance of CsH" 1 ". The illuminating star is positioned at the right of the plots. The symbols 
present the measured range of possible abundances for C2H (blue vertical segment), C3H2 (green vertical segment), and 
CaH + (red filled rectangle), inferred at a typical resolution of 6". The legend at the bottom of the figure presents the 
line coding for the curves: C2H in plain blue, C2H" 1 " in dashed blue, C3H in plain red, C3H" 1 " in dashed red, and C3H2 in 
green. 



We used the Ohio State University (osu) pure gas- 
phase chemical network upgraded for photochemical stud- 
ies. CsH + is produced by reactions between C2H2 and 
C+. Then, C 3 H+ is thought to produce C 3 H+, and C 3 H+ 
through reactions with H2, which later recombines with 
electrons to form C3H and C3H2 

h 2 c 3 H+^-^C 3 H 

c + 4. - 

C2H2 — 5- C3H 

^ C 3 H+ -U- C 3 H 2 . 

Hence these species are usually included in gas-phase chem- 
ical networks with reaction rate accuracies of a factor 2 or 
better. We can thus use them to compare models with ob- 
servations. 

FigureOpresents the results of the photochemical model 
for a few small hydrocarbon molecules and cations, namely 
C 2 H, C 3 H, C3H2, C 2 H+, and C 3 H+. From top to bot- 
tom, the figure shows the spatial profiles of the density, 
abundances relative to the number of hydrogen atoms, and 



the abundances relative to CsH + . The left column presents 
the profiles computed by the code, which samples the UV- 
illuminated gas on a finer spatial grid than the UV-shielded 
gas to correctly represent the steep physical and chemical 
gradients. The right column presents the profiles convolved 
with a Gaussian of 6" full width at half maximum to fa- 
cilitate the comparison with the abundances inferred from 
PdBI observations at 6" angular resolution. The measured 
abundances are displayed with vertical segments for C2H, 
and C-C3H2, and with a filled rectangle for C3H" 1 ". Although 
the 1-C3H+ abundance is only inferred from the IRAM-30m 
single-dish telescope, we also show it here because we used 
a simple model of its emission based on the PdBI hydrocar- 
bon data to correct for the beam dilution. Table [3] quanti- 
tatively compares the measured and modeled abundances 
at the PDR position ("IR peak" at A v ~ I mag) , inside 
the molecular cloud at a position representative of the UV- 
shielded material "cloud" ( A v ~ 8 mag) , and closer to the 
Hn region (i.e., "IR edge" at A v ~ 0.01 mag). 
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We only consider the most stable isomers to compare 
the measured and modeled abundances. Table [4] indicates 
that the amount of C3H2 locked in the linear species is 
negligible, and that twice as much C3H is locked in the 
cyclic spe cies compared to the linear species. Moreover, 
Table 1 of lMebel et al.l (|2007l) indicates that the cyclic form 
of C3H2 and C3H are more stable than their linear form 
by 59.3 Llmo!" 1 (or 0.61 eV), and 11 kJ mol" 1 (or 0.11 eV), 
respectively. Finally, ISavic et all (|2005l ) indicated that the 
cyclic form of CsH + is less stable than the linear one by 
220kJmol~ 1 , or 2.28 eV. Given the large difference in en- 
ergies, we here assume that the linear form is the main 
product in the gas phase. 

While C2H has much brighter lines and a column den- 
sity higher by two orders of magnitude compared to linear 
and cyclic C3H, we have no hint of a detection of bright 
lines associated to C2H + . This is due to a combination 
of two effects. First, the spectroscopic structure of C2H+ 
(which has a 3 I1 ground electronic state) is more complex 
than that of C3H + , implying that the emission is spread 
over more lines. The knowledge of the rotational spectrum 
of C2H + would enable a deep search of this reactive molec- 
ular ion to check whether C2H" 1 " is also present in the horse- 
head PDR. Second, C2H" 1 " is predicted to be less abundant 
than C3H + (except when A v < 0.1 mag), because it is effi- 
ciently destroyed by H2 , independently of the temperature. 
In contrast, the detection of C3H+ at the warm PDR po- 
sition corroborates that the destruction of C3H + by H2 is 
much less efficient there than in the cold-core region, i.e., 
the C3H + + H2 reaction is strongly dependent on the gas 
temperature ([Savic fc G erlich 20051). 

On one hand, the agreement between the interferomet- 
ric abundances of C2H and C3H2 and the modeled curves 
improved at the "cloud" position by up to four order s of 
magnitude compared to the models of lPetv et al.l (|2005l . see 
their Fig. 10). This is the result of the recent addition of the 
chemistry on the grain surface to the Meudon PDR code. 
Indeed, the modeled and predicted abundances of C-C3H2 
are consistent, while the modeled abundance of C2H is still 
one order of magnitude lower than the measured abundance 
(cf. Table [2]). However, it is worth noting that the "cloud" 
position (Sx — 24.9", Sy — —5.3") does not fall, for histor- 
ical reasons, on the dense core. Indeed, the C2H emission 
presents a dip at the DCO+ peak (Sx = 44.7", Sy = 16.5", 
blue cross on Fig. [2]), while the C-C3H2 emission shows a 
bright filament there. These observational facts indicate 
that C2H probably depletes at a higher rate onto grains 
than C3H2. 

On the other hand, the model results did not change 
significantly in the UV-illuminated part (IR peak and IR 
edge). The measured abundances are typically one and 
two orders of magnitude higher than the modeled abun- 
dances for C2H and C3H2. In contrast, the measured abun- 
dance of C3H+ is only about a factor 2 higher than the 
modeled abundance (cf. Table [3]) . If the assumptions used 
here are correct (the geometry of the I-C3H 4 " emission, 
and neglecting the amount of c-CsH + ), this difference be- 
tween the small hydrocarbon molecules and 1-C3H + would 
confirm that there is a non gas-phase chemical way in 
which small hydrocarbons are formed i n UV-illuminate d 
regions, e.g., the photo-erosion of PAHs dPetv et al.ll2005lh 
Interferometric observations of C3H and C3H" 1 " are needed 
to confirm this result. Indeed, Fig. [5]shows the spatial pro- 
file of the modeled C3H 4 " abundance convolved by convo- 
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Fig. 6. Spatial profiles of the predicted abundance of CsH + 
cation. The modeled abundance profile was convolved along 
the x-axis with a Gaussian of 6"-FWHM (red line), and 
27"-FWHM (blue line). 

lution at 6" and 27". Clearly, only interferometric obser- 
vations can provide the high angular resolution needed to 
resolve the actual structure of the C3H 4 " emission, and po- 
tentially detect a spatial shift with respect to the millimeter 
wave emission of other small hydrocarbons and to the in- 
frared imaging of PAHs and dust grains. 

6. Summary 

We reported the first detection of 1-C3H+ in the interstel- 
lar medium. Laboratory measurements of the 1-C3H+ spec- 
troscopy are needed to improve the spectroscopic character- 
ization of this molecular ion. Interferometric imaging with 
either PdBI or ALMA is required to better constrain the 
small hydrocarbon chemistry. 
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Appendix A: Fit likelihood 

We wish to know to which development order we can sig- 
nificantly fit the frequency model to the observed data. 
Assuming that all intensity channels are independent mea- 
sures with Gaussian noise, the likelihood, £, of the fit of 
the model to the data is given by 



where i is the index over the N measured channel intensi- 
ties, Oi the noise of channel i, I° hs and 7 mod the observed 
and modeled intensities, and 9 the set of model parameters. 
Taking the logarithm yields ln(£) = K — 0.5x 2 , where 



K is a constant and x 2 is the usual chi square measure of 
the fit quality. If we try to fit two models M\ and M 2 , 
which differs only by one parameter, to the same data, 
the comparison of the fits is said to have only one degree 
of freedom. The difference x'm 1 ~ Xm 2 = 2 In (Cm 2 I ) 
then follows the well known x 2 distribution with one de- 
gree of freedom. In our case, the difference of % 2 between 
the second and third order models is 9.2, which implies that 
the third-order model better reproduces the data than the 
second-order model with a probability > 99.75%. In con- 
trast, developing the centrifugal effect to the fourth order 
does not bring a significant improvement over the third or- 
der development. We therefore stopped the determination 
of the spectroscopic constants to the third order model. 




(A.l) 




N 
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Fig. A.l. Hydrocarbon lines detected toward the Horsehead PDR at RA=5h40m53.936s, Dec=-2°28'00" (J2000). Only 
lines in the 3 mm band, whose energy of the upper level is lower than 100 K and which have a higher signal-to-noise ratio 
than 8, are shown. The green lines are Gaussian fits. The vertical dashed lines mark 10.7 fans -1 . 



Table A.l. Gaussian fit results for the lines of small hydrocarbons detected toward the Horsehead PDR at 
RA=5h40m53.936s, Dec=— 2°28'00" (J2000). Only lines in the 3mm band with an energy of the upper level lower 
than 100 K and a higher signal-to-noise ratio than 8 are shown. 



Species 






Transition 






Frequency 
MHz 


Area 
mK km s~ 


Velocity 
kms -1 


Width 

kms -1 




N 


= 1 


-0, J = 3/2- 1/2 F = 1- 


1 




87284.105 


238±3 


10.675±0.005 


0.735±0.013 




N 


= 1 


-0,J = 3/2-l/2F = 2- 


1 




87316.898 


1943±6 


10.682±0.001 


0.815±0.003 


C 2 H 


N 


= 1 


-0, J = 3/2 -1/2 F = 1- 







87328.585 


1035±4 


10.675±0.002 


0.769±0.004 


N 


= 1 


-0,J = 1/2-1/2 F = 1- 


1 




87401.989 


1103±5 


10.694±0.002 


0.785±0.004 




N 


= 1 


-D,J = l/2-l/2F = 0- 


1 




87407.165 


489±4 


10.697±0.003 


0.749±0.008 




N 


= 1 


-0, J= 1/2- 1/2 F= 1- 







87446.470 


262±3 


10.700±0.004 


0.745±0.011 




N = 


2l2 


- In, J = 5/2 - 3/2 F = 3 




2 


91494.231 


106±3 


10.395±0.010 


0.691±0.024 


c-C 3 H 


N = 


2l2 


- In, J = 5/2 -3/2 F = 2 




1 


91497.525 


59±3 


10.561±0.012 


0.594±0.034 


N = 


2l2 


- In, J = 3/2 - 1/2 F = 1 







91692.752 


32±3 


10.712±0.036 


0.887±0.084 




N = 


2l2 


- In, J = 3/2 - 1/2 F = 2 




1 


91699.471 


74±4 


10.760±0.020 


0.858±0.050 






J = 


= 9/2-7/2 F = 5-4(/) 






97995.166 


44±7 


10.605±0.051 


0.647±0.108 


1-CsH 




J: 


= 9/2 - 7/2 F = 4 - 3(/) 






97995.913 


45±10 


10.583±0.092 


0.848±0.235 




J 


= 9/2 - 7/2 F = 5 - 4(e) 






98011.611 


41±7 


10.629±0.046 


0.599±0.116 






J 


= 9/2-7/2 F = 4-3(e) 






98012.524 


35±8 


10.550±0.075 


0.651±0.183 








2q 2 — 111 






82093.542 


242±6 


10.704±0.009 


0.741±0.023 








3l2 — 3o3 






82966.200 


68±3 


10.660±0.012 


0.618±0.029 


C-C3H2 






322 — 3l3 
2l2 — loi 
3o 3 — 2i2 






84727.696 
85338.893 
117151.191 


14±2 
748±4 
507±72 


10.633±0.028 
10.687±0.002 
10.680±0.051 


0.344±0.051 
0.767±0.005 
0.715±0.117 








4l4 - 3l 3 






82395.090 


107±7 


10.681±0.020 


0.720±0.057 








4o4 — 3o3 






83165.256 


43±4 


10.362±0.029 


0.647±0.068 








4l 3 — 3l2 






83933.700 


93±3 


10.690±0.010 


0.689±0.025 


1-C 3 H 2 






5l5 — 4i 4 

5 33 — 432 

505 — 4o4 
5l4 — 4i3 






102992.379 
103914.354 
103952.926 
104915.583 


81±3 
27±2 
34±3 
58±4 


10.658±0.012 
10.592±0.023 
10.803±0.030 
10.696±0.025 


0.657±0.030 
0.538±0.052 
0.652±0.068 
0.703±0.058 






N 


= 9-8,J = 19/2-17/2 






85634.004 


198±4 


10.737±0.007 


0.689±0.018 






N 


= 9-8, J = 17/2-15/2 






85672.579 


168±4 


10.725±0.007 


0.663±0.018 






N 


= 10-9, J = 21/2 -19/2 






95150.388 


204±3 


10.734±0.004 


0.682±0.011 


C 4 H 




N 


= 10 -9, J = 19/2 - 17/2 






95188.946 


177±3 


10.700±0.005 


0.665±0.013 




N = 


= 11 - 10, J = 23/2 - 21/2 






104666.565 


202±4 


10.732±0.007 


0.714±0.019 






N = 


= 11 - 10, J = 21/2 - 19/2 






104705.106 


177±5 


10.705±0.011 


0.728±0.026 






TV = 


= 12 - 11, J = 25/2 - 23/2 






114182.512 


162±12 


10.665±0.021 


0.582±0.051 






N = 


= 12- 11, J= 23/2 - 21/2 






114221.040 


186±12 


10.663±0.022 


0.719±0.056 



